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We study numerically the cooling of a young bare strange star and show that its thermal lumi- 
nosity, mostly due to e'''e~ pair production from the quark surface, may be much higher than the 
Eddington limit. The mean energy of photons far from the strange star is ~ 10^ keV or even more. 
This differs both qualitatively and quantitatively from the thermal emission from neutron stars and 
provides a definite observational signature for bare strange stars. It is shown that the energy gap 
of superconducting quark matter may be estimated from the light curves if it is in the range from 
~ 0.5 MeV to a few MeV. 
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Strange stars that are entirely made of deconfined 
quarks have been long ago proposed as an alternative to 
neutron stars (e.g., jl)). The bulk properties (size, mo- 
ment of inertia, etc.) of strange and neutron stars in the 
observed mass range (1 < M/Mq < 2) are rather simi- 
lar, and it is very difficult to discriminate between strange 
and neutron stars . Strange quark matter (SQM) with 
a density of ~ 5 x 10^"' g cm^'^ might exist up to the sur- 
face of a strange star [Q. This differs qualitatively from 
the neutron star surface and opens observational possi- 
bilities to distinguish strange stars from neutron stars. 

Normal matter (ions and electrons) with a mass AM < 
10~^Mq may be present at the SQM surface of a strange 
star 1^ . Such a massive envelope of normal matter could 
completely obscure the quark surface. However, a strange 
star at the moment of its formation is very hot. The 
temperature in the stellar interior may be as high as a 
few times 10^^ K |^. The rate of neutrino-induced 
mass ejection from such a hot compact star is very high 
1^. Therefore, in a few seconds after the star formation 
the normal-matter envelope is blown away, and the SQM 
surface is nearly (or completely) bare 0. A strange star 
remains nearly bare as long as the surface temperature 
is higher than - 3 x lO"^ K g. 

Due to the high plasma frequency of SQM, Up ~ 20 
MeV, a bare strange star will be a very inefficient emit- 
ter of thermal X-rays as soon as its temperature drops 
below ^ 10^^ K, i.e., a few seconds after it is born [^. 
This fact suggested that bare strange stars would be very 
difficult to detect. However, the enormous surface elec- 
tric field, which binds the electrons to the quark matter, 
will induce intense emission of e"'"e~ pairs |^] and sub- 
sequent hard X-ray emission, at luminosities above the 
Eddington limit, LEdd ^ 1-3 x 10^^{M/Mq) erg s~\ as 
long as the surface temperature is above ^ 5 x 10^ K 
0. This process significantly increases the possibifity 
to detect young bare strange stars, but its importance 
depends on how long a high luminosity can be sustained. 

We want to address this issue in this letter by modeling 
in detail the thermal evolution of a young bare strange 
star in order to calculate its light curve, considering var- 



ious scenarios about the state of SQM. It should be em- 
phasized that the resulting evolution described here will 
differ both qualitatively and quantitatively from the evo- 
lution of a strange star with a crust, i.e., covered by a 
small layer of normal matter, [ p^ , and from the evolution 
of a more standard compact object ]ri[ . 

Recently, it has been argued that SQM is a color 
superconductor with a very high critical temperature 
Tc ~ 10^^ K (for reviews, see |l^). At extremely high 
density the color superconductor is in the "Color-Flavor- 
Locked" (CFL) phase in which quarks of all three fiavors 
and three colors are paired in a single condensate. In 
this CFL phase SQM is electrically neutral and no elec- 
trons are present. If the strange quark mass, m^, is not 
too large the CFL phase may extend down to the lowest 
density corresponding to the surface of a strange star, in 
which case the considerations of this paper are irrelevant 
since no electrons would be present at the surface and 
hence there would be no supercritical electric field and 
no e"'"e~ pair emission. However, for sufficiently large rus 
the low density regime is rather expected to be in the "2 
color-flavor SuperConductor" (2SC) phase in which only 
u and d quarks of two color (say ui, U2, di and ^2) are 
paired in a single condensate while the ones of the third 
color, say U3 and c?3, and the s quarks of all three colors 
are unpaired. In the 2SC phase electrons are present. We 
will consider strange stars made entirely of normal, i.e., 
unpaired, SQM, which may be unrealistic but is a bench- 
mark and the case of SQM in the 2SC phase only and 
with a mixture of 2SC phase with CFL phase at high den- 
sity and, finally, SQM in the 2SC phase with secondary 
pairing of the u^-d^ and s quarks with a much smaller 
gap. In all cases the critical temperature is related to the 
energy gap A(0) at T = through ~ 0.57A(0) as in 
BCS theory (in natural units h = c = ks = 1). 

We consider, as a typical 1.4 Mq strange star 

which we construct by solving the Tolman-Oppenheimer- 
Volkoff equation of hydrostatic equilibrium using an 
equation of state for SQM as described in |j [with 
a bag constant B = (140MeV)'*, QCD coupling constant 
as = .g^/47r = 0.3, and rUg = 150 MeV]. The thermal 
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evolution is determined by the energy conservation and 
heat transport equations: 
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where is the specific heat of the matter, n its thermal 
conductivity and Qi, its neutrino emissivity, and Fr is the 
heat flux at radius r. We are actually solving the general 
relativistic version of these equations with a Henyey-type 
cooling code jl^. At the stellar surface, the heat flux di- 
rected outward the strange star is equal to the energy flux 
emitted from the stellar surface, Fr{r ^ R) = F-^ + F±, 
where F^ and F± are the energy flux emitted from the 
unit surface of SQM in thermal photons and e^e~ pairs, 
respectively. This and equation give the boundary 
condition on dT/dr at the stellar surface. The bound- 
ary condition at the stellar center is dT/dr = 0. In our 
numerical simulations, we adopt the values of F-y and 
F± from 0. We assume that at the initial moment, 
t = 0, the temperature in the strange star is uniform, 
T = 10" K - 10 MeV. 

The specific heat is Cv — Cy{e) + J2q^vil) (? 
ning through the nine quark components, ui, S3) with 
Cv{i) = ^IJ'iT, fii being the chemical potential of the 
ith component. When a component i becomes paired 
its Cvii) first increases by a factor 2.426 at Tc and then 
decreases exponentially at T -C Tc Q. 

The neutrino emission is due to the three quark direct 
Urea processes Uc + e~ — > dc + Ve and dc — > Wc + e~ -t-Fe, 
c = 1, 2, 3 denoting the quark color which is not altered 
by weak interactions, with emissivities [151 
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and the corresponding, but weaker, processes with dc ^ 
Sc- These processes are also strongly suppressed when 
one of the participating component is paired . 

The thermal conductivity is an essential ingredient of 
our calculations since it determines how fast heat is car- 
ried to the surface from the underlying layers. When no 
quark pairing is present k has been calculated in mTj : 
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In the case of the 2SC phase the thermal conductivity 
will be provided dominantly by the unpaired quarks 1*3 
and which only suffer scattering through the exchange 
of the fully screened (both Debye and Meissner) gluon of 
adjoint color index 8 [^. Following the method of ||l^ 
we obtain 
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In the CFL phase the thermal conductivity is extremely 
large [p^ because this phase is transparent to photons, 
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FIG. 1: Upper panel: luminosities in neutrinos, L^, and pho- 
tons plus e'*'e~ pairs, Lth, for a strange star in the normal 
phase (solid lines) and in the 2SC phase (dashed lines). Lower 
panel: temperature in the center, Tc, at the surface, T, and 
1 meter below the surface, Ts-im of the strange star in the 
normal phase. The large difference between T and Ts-im 
at early times illustrates the enormous temperature gradient 
present just below the surface because of the large Lth. 



Figure |^ shows the temperature and luminosities (in 
both neutrinos and surface thermal emission) as a func- 
tion of time t, for a strange star in the normal phase 
and in the 2SC phase with a gap A2SC ^ 100 MeV. The 
neutrino luminosity, L^, is always much higher than the 
surface thermal luminosity, ith, i-e., neutrinos drive the 
cooling and the surface emission just follows the evolu- 
tion of the bulk of the star. The occurrence of the 2SC 
phase has little effect on the star's evolution: the rea- 
son is that the neutrino emission is cut by a factor three 
since u and d quarks of color 1 and 2 do not contribute 
anymore, and the specific heat is cut by a factor 9/5. 
The neutrino cooling time scale is TcooI ~ CyT/Qi, (see 
Eq. (|l|)), and it therefore does not change significantly. 

We see that the thermal luminosity may be many or- 
ders of magnitude higher than the Eddington limit for a 
period of several days after the strange star formation. 
Such a high luminosity is allowed for a bare strange star 
because at its surface SQM is bound via strong interac- 
tion rather than gravity and, therefore, a bare strange 
star is not subject to the Eddington limit in contrast to 
a neutron star ^. Super-Eddington luminosities are 
a fingerprint of hot bare strange stars. At t > 10^ s, 
the thermal emission decreases fast, and after about one 
month when < 3 x 10^ K, the thermal radiation be- 
comes practically undetectable (see Fig. |^). This strong 
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decrease of the surface emission is due to the suppression 
of the e"''e~ pairs emission by increasing degeneracy of 
electrons in the thin "electron atmosphere" 0. 

Assumption that a part of the interior of the strange 
star is paired in the CFL phase, with Acfl > 10 MeV, 
has no effect at all on the surface thermal luminosity. The 
reason is simply that it cuts down both and Cy in the 
same fraction and does not affect TcooI at all. We have 
checked it numerically with a mass of up to 1.39 M0 of 
the star, for a total mass of 1.4 Mq , in the CFL phase 
and found variations in Lth smaller than the thickness of 
the line in the plot of Fig. |l|. The star's temperature pro- 
file is of course affected since the CFL core cools mostly 
by heat diffusion into the outer 2SC phase layer, but neu- 
trino emission in this outer layer is so strong that Tg and 
Lth are practically not affected. This result is similar to 
the case of quark matter in neutron star interiors which 
is undetectable if it is in the CFL phase 0. 

Two effects can dramatically alter the results of Fig. |^. 
They are a secondary pairing of the M3 and quarks, 
and also possibly of the s quarks, in the 2SC phase, and 
convection in the upper layers |Q . 

Pairing of the M3 and quarks, in addition to the 2SC 
pairing, will quench the 3rd color channel of neutrino 
emission by the direct Urea process of Eq. ^ and reduce 
the specific heat while pairing of the s quarks will have 
little effect on the neutrino emission but will reduce C„. 
The suppression of Qi, and Cy is very strong only if the 
resulting gap is nodeless, while a gap with nodes on the 
Fermi surface will only produces a reduction of the order 
of {T/TcY if nodes are at isolated points and of order of 
T/Tc if nodes are ID lines. The coohng time scale, TcooI, 
can be increased or reduced depending on whether Qi, 
or Ctj is the most strongly suppressed. Since very little 
is known on these possible secondary gaps (see however 
pl[) we consider three cases for the 1*3-^3 gap A3 and 
the s gap A,.: [A] A3 = 0.3 Mev and A^ = 3 MeV, [B] 
A3 = As = 3 MeV, and [C] A3 = 3.0 Mev and A,. = 0.3 
MeV. These three cases span a large range of suppression 
of Qn and/or 

Figure || shows the resulting range of Lth possibly 
attained through this secondary pairing. We assumed 
nodeless gaps to maximize the effects, but have verified 
explicitly that gaps with nodes give results which are in- 
termediates between the three cases presented here. Case 
[A] has very little Qn suppression and only partial reduc- 
tion of Cv , and it naturally differs little from the case of 
no pairing. Cases [B] and [C] have strong suppres- 
sion and their cooling is eventually driven by the surface 
thermal emission. Since case [B] also has a strong C„ 
reduction it cools faster than case [C] whose Cy is more 
moderately suppressed. The knees on the curves [B] and 
[C] at times t ^ 10'' sec and t ^ 3x 10^° sec, respectively, 
correspond to the moment when the star becomes almost 
isothermal. 

In all the models presented above, at early times, the 
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FIG. 2: The thermal luminosity of the strange star in the 2SC 
phase with secondary pairing of the Ms-da and s quarks for 
three different scenarios of pairing [A] , [B] and [C] as described 
in the text (solid lines). Arrows on curves [B] and [C] point 
the time at which L^, becomes lower than Lth', in case [A] 
ii/ ^ Lth at all times. The dotted line shows the evolution 
of Lth without any pairing for comparison (from Fig. |l|). 



temperature gradient just below the surface (see Fig. ||, 
lower panel) is many orders of magnitude higher than the 
adiabatic temperature gradient p2[ 
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where riq is the quark number density. If convection can 
develop, given the shallowness of the superadiabatic layer 
we can, as a first approximation, consider the star as hav- 
ing a uniform temperature. Moreover, in the case of pair- 
ing of the M3-c?3 quarks, and/or s quarks, we may expect 
convective counterflow of the u^-ds, and/or s, superfluid 
which may be even more efficient that convection to erase 
any temperature gradient. 

Figure || shows the resulting range of Lth for isother- 
mal stars for the same three cases of secondary pairing 
as in Fig. |^. In cases [B] and [C], a few seconds after 
the star formation the u^-d^ pairing occurs and subse- 
quently L^ < Lth while in case [A] neutrino losses drive 
the cooling during all times. The sharp drop of Lth in 
case [B] at t ~ 20 s is due to the strong suppression of 
Cv from pairing and does not occur in cases [A] and [C] 
in which the u^-d^ or s gaps are smaller. At early times 
the isothermal models have a much higher thermal lumi- 
nosity that the diffusive ones of Fig. |^ since their surface 
temperature is much higher. However, the isothermal 
models whose cooling becomes driven by thermal emis- 
sion, [B] and [C], cool faster and eventually have a lower 
Lth than the diffusive ones. Naturally, isothermal and 
diffusive models follows the same evolution once the lat- 
ter become isothermal, i.e., after the knees of cases [B] 
and [C] in Fig. |. 

How much the density dependence of the chemical 
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ton spectrum together with the hght curves calculated in 
this paper could be a good observational signature of a 
young bare strange star. 
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FIG. 3: The thermal luminosity for the same scenarios as in 
Fig. ^. It is assumed that the star is isothermal due to either 
convection or superfluid counterflow. The dotted curve shows 
the evolution of Lth without any pairing for comparison. 



composition of SQM, and of the pairing gaps, may re- 
duce convection, and superfluid counterflow, is an open 
question and precludes us to firmly decide which cooling 
trajectories, Fig. || or ^, are the correct ones. However, 
the gradients of both /i^'s and A's are small in SQM and 
one may expect only small reductions, i.e., the isother- 
mal models of Fig. ^ are probably more appropriate that 
the diffusive ones of Fig. ||. 

In our simulations we assumed that neutrinos escape 
freely from the stellar interior. This is valid only in a 
few seconds after the strange star formation when the 
internal temperature is less than ~ 10^° K dl- In this 
case, e+e" pairs created at the surface of SQM prevail 
over photons in the surface thermal emission H but pairs 
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pho lons in the viciuiL^' uf Lhe sLiaugu star 83 

In the process of the star cooling the photon spec- 
trum varies significantly. At very high luminosities, 
Lth > 10''^ ergs s~^, we expect that the photon spec- 
trum is nearly blackbody with a temperature Tbb — 
To(ith/10'*3 erg s-i)!/^^ where Tq - 2 x 10^ K Q. For 
intermediate luminosities, 10''^ < Lth < 10'*'^ erg s^^, the 
effective temperature of photons is more or less constant, 
Tbb To 11. At Lth < 10"*^ erg s~i, the photon spec- 
trum essentially differs from the blackbody spectrum, 
and its hardness increases when Lth decreases. This is be- 
cause photons that form in annihilation of e+e^ pairs do 
not have enough time for thermalization before they es- 
cape from the strange star vicinity. When the photon lu- 
minosity decreases from ~ 10*^ erg s^^ to ^ 10'^^ erg s^^, 
the mean energy of photons increases from ~ 30 keV 
to ~ 500 keV while the spectrum of photons eventually 
changes into a very wide {AE/E ^ 0.3) annihilation line 
of energy E 500 keV . Such a variability of the pho- 
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